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Correlating features in the primordial spectra 
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Heavy fields coupled to the inflaton reduce the speed of sound in the effective theory of the 
adiabatic mode each time the background inflationary trajectory deviates from a geodesic. This can 
result in features in the primordial spectra. We compute the corresponding bispectrum and show 
that if a varying speed of sound induces features in the power spectrum, the change in the bispectrum 
is given by a simple formula involving the change in the power spectrum and its derivatives. In this 
manner, we provide a uniquely discriminable signature of a varying sound speed for the adiabatic 
mode during inflation that indicates the influence of heavy fields. We find that features in the 
bispectrum peak in the equilateral limit and, in particular, in the squeezed limit we find considerable 
enhancement entirely consistent with the single field consistency relation. From the perspective of 
the underlying effective theory, our results can be generalized to incorporate a wide variety of 
inflationary models where features are sourced by the time variation of background quantities. 



Although no significant evidence for features in the pri- 
mordial power spectrum has been observed to date [I], 
a 10% modulation of power remains entirely consistent 
with direct reconstruction [5] of that region of the power 
spectrum accessible to present cosmic microwave back- 
ground (CMB) observations [3 . Further improvement in 
our knowledge of the matter spectrum beyond the CMB 
scales from future large scale structure surveys [4j ought 
to furnish far superior statistics that may yet allow us 
to determine if there are any other scales present in the 
primordial power spectrum other than that which sets its 
amplitude [5]. 

Inflationary models that generate features have been 
considered in the literature for various theoretical moti- 
vations. Examples include explaining outliers to the best 
fit ACDM models in the CMB [5] (see however [7]), prob- 
ing and constraining couplings of the inflaton to other 
fields through particle production [8], probing modified 
vacuum structure [9], interrupted slow-roll |10) and per- 
haps even signatures of inflation's embedding in a UV 
complete theory [TT] . 

The perspective of this article is informed by recent 
findings that features induced by reduced and varying 
speeds of sound c s of the adiabatic mode might be an 
inevitable consequence of field excursions during infla- 
tion [12HT7] . and the surprisingly large role played by 
heavy fields in the dynamics of light fields (first noticed 
in [T5], additionally elaborated upon in [IH])- Specifically, 
heavy fields with masses much larger than the Hubble 
scale can influence the dynamics of the adiabatic mode 
in a manner that is entirely consistent with the persis- 
tence of slow- roll [HUE], the validity of an effectively sin- 
gle field description and the decoupling of the true high 
and low energy modes of the system [16] , which further- 



more remains weakly coupled up to the cutoff defining 
the theory 20J . At low energies, the resulting effective 
theory of the adiabatic mode is an operator expansion 
parametrized entirely by c s [15], which can in fact be a 
relatively rapidly varying function, subject to the adia- 
baticity condition discussed below. 

Here we show that the generic consequence of a vary- 
ing c s would be to imprint features in the power spec- 
trum and enhance the bispectrum for all shapes, peak- 
ing in the equilateral but also enhancing other config- 
urations (see |21j for similar studies). Our main result 
is summarized by|l3L where / NL parametrizes the bis- 
pectrum as in ( |10[ ), in a particular configuration A of 
the three wavevectors. There, the non-trivial scale de- 
pendence of is directly correlated to the features in 
the power spectrum parametrized by AT-n/Tii, while the 
c A -coefficients depend only on the shape of the specific 
configuration A. This enhancement of the primordial 
bispectrum in combination with correlated features in 
the power spectrum could be readily discernible through 
future observations. 

Our starting point is the effective action for the adi- 
abatic mode. The quadratic action S2 obtained by in- 
tegrating out a heavy field in a two-field model is given 
by 



S 2 = 



'pi 
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where a is the scale factor, H = a/a and e = —H/H 2 . 
Here, c s is the speed of sound of the adiabatic fluctuations 
given by [15] 
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where M 2 S = M 2 — 9 2 is the effective mass of the nor- 
mal (isocurvature) fluctuations off the background tra- 
jectory, is the angular velocity of the background in 
field space [T3] , and M 2 is its mass 2 in the absence of 
any turns. A particular situation of interest is when the 
background trajectory departs off the adiabatic minimum 
of the potential as inflation progresses due to 'bends' in 
field space, resulting in transient reductions in the speed 
of sound [25] HH 113] • Presuming that this reduction in 
the speed of sound is sufficiently small (|1 — c~ 2 | <C f) 
we conveniently rewrite the quadratic part of the action 
^ as 



S2 — SWree + 



74 S 2 

a xa enipi 



ii 2 



(3) 



where S^froc is nothing but S2 with c s — 1. With this 
splitting of the quadratic action, S^.mt may be consid- 
ered as a perturbation to $2, free- Thus one can readily 



evaluate the corrections to the power spectrum and the 
bispectrum induced by a varying c s . Denoting |22| 
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the change in the power spectrum generated by changes 
in the speed of sound, to first order in u, is given by [53] 
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dr u{t) sin (2fcr) 
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where V-jz = H 2 / (8ir 2 em 2 ,^ is the featureless flat power 
spectrum and r is the conformal time. One immediately 
sees how features in the power spectrum are generated 
by a varying c s . We also note that constant, reduced 
speeds of sound renormalize the power spectrum without 
imparting any new scale dependence. 

One can also compute the leading contributions to the 
bispectrum AB-r due to a non- vanishing u, as |23j 



AB n {k x , k 2 , fc 3 ) = 23* <^ 2^(0)^(0)^(0) 
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where the mode function solution IZk (r) is given by 



n k ( T ) = 



iH 



%AhF 



m P1 



(1 + ikr) e 



— ik7 



(7) 



In the above, we work to leading order in slow-roll param- 
eters and to linear order in u. The complete expression 
for the bispectrum is of the form B = Bq + AB, where 
Bq represents the leading non-zero contributions when 
c s = 1 and is of 0(e 2 ). AB dominates whenever u's 
maximum value |it| max is larger than 0(e). From ^ we 
see that |M| max of O(10 _1 ) translates into features in the 
power spectrum of 0(10)%, which is reasonable for the 
scales accessible to the CMB to this level. Thus, with 
slow-roll parameters no bigger than 0(1O~ 2 ), we obtain 
that AB would become the dominant contribution to the 



bispectrum. 

Taken separately, one immediately sees from ([5| and 
^ how the changes in c s parametrized by u source scale 
dependence in both the power spectrum and the bispec- 
trum. However to highlight their correlation in a man- 
ner that should serve as a useful discriminant in probing 
data, it is useful to invert u to linear order in terms of 
the change in the power spectrum as 



u{t) = 
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where u is defined as the odd extension of u over the 
real line. Substituting this into ^ allows us to infer the 
leading contribution to the bispectrum as 
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where = |fci| and k = (ki + k 2 + k 3 )/2 [23] ■ This is the result we wish to highlight: features in the power spectrum 
translate directly into correlated features in the bispectrum, with the precise nature of the correlation depending on 
the configuration we look at. We may define a dimensionless shape function, with the /nl ansatz in mind, as 
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where the shape superscript on the left hand side indicates that it is a function of particular configurations of the 
three wavevectors. Evaluating the above in certain interesting limits, we find 
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We notice that the squeezed limit is nothing other than a 
re-expression of the single field consistency relation [2"4"] 



Bn(h,h,h) 



dlogPn(k 3 ) 
d log k 3 
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with P n {k) = 2TT 2 Vn/k 3 . Thus may approach 

values of order unity, even though slow-roll is operative 
throughout. This is because the spectral index receives 
contributions that go like s = c s /(Hc s ) [5S] which can 
approach order unity consistent with slow-roll and the 
validity of the single field approximation [T3TU6| . From 
(111 we infer that in general, for any configuration one 



might look at, the bispectrum is sourced by features in 
the power spectrum as 



/nl ~ c o (k) 
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with a prime denoting a logarithmic derivative and the 
coefficients cf" depending only on the shape of the con- 
figuration we look at, and with all information about 
the variation of c s encoded in AVn/Vn- In Figure [l] 
we plot the changes in the power spectrum against the 
bispectrum for a few prototypical changes in the speed 
of sound that model different varieties of turns in the 
background trajectory. As evident from the plots, the 



shape function in each limit ( 11 ) exhibits features in ac 



cordance with modulations in the power spectrum and 



its derivatives with the amplitude of the equilateral /i- eq ' ) 



'NL 



being as high as ~ 3 for the parameters we consider if we 
take | | max ~ 10 _1 . More interestingly, the squeezed 
configuration, in addition to being ir/2 out of phase 
with the modulations of the power spectrum as expected 



f( s q) 



from the consistency relation (12 1, attains a peak value 
0.5. Although this is very unlikely to be de- 



tectable in the CMB, the prospects for measuring /nl to 
a precision of 0(1) at smaller scales stands to improve 
with future observations [4]. 

From the perspective of the effective field theory of 
inflation of [26] , the cubic order action is parametrized 
by coefficients M 2 and M 3 in such a way that c s is 
determined by M 2 . On general grounds [27] one ex- 
pects M 3 to be of the form M| - (l - c~ 2 ) M|, with 
the precise relation encoding properties of the parent 
theory from which the effective theory descends. For 



'■) for DBI inflation, 



example, M|/M 2 4 = -3(1 - 
whereas M 3 4 /M 2 4 = 3(1 - c7 2 )/4 when one heavy field 
has been integrated out in a two-field theory [15]. Be- 
cause M 3 ~ (l — cj 2 ) M 2 implies terms of 0(u 2 ) in |6l, 
we see that ^ and ( [TT| ) are representative of a wide 
class of single field models of inflation, beyond the spe- 
cific case in which adiabatic modes interact with heavy 
fields. In spite of these details, our results show that one 



will always arrive at the general expression ( 13 ) where 
the coefficients c~~ depend on the details of the parent 
theory as well as the shape of the configuration we are 
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FIG. 1: AVm/Vn (solid line) versus (equilateral, dotted line, left panel) and /^iP (squeezed, dotted line, right panel) 

for u = -u max {tanh[(iV-7Vi)/<WV] - tanh \(N - N f )/SN}} with u max = l/12,N f - Ni = 2 and 5N = 0.05 (top panel), and 
u = — u max / cosh 4 [(N — Ni,)/AN] with u max = 1/12 and AN — 0.2. TV = — log|r| is the number of e-folds, ★ denotes a 
convenient reference, and here the central moment of the turn, and and k+ are the corresponding values. 



interested in. 

The prospects for observing a non-trivial scale depen- 
dence in the primordial bispectrum are enhanced if they 
are correlated at commensurate comoving scales with fea- 
tures in the power spectrum. In this article, we illustrate 
a context in which this occurs naturally- when the dy- 
namics of the adiabatic mode are influenced by heavy 
fields in such a way as to transiently reduce c s at var- 
ious points along the inflationary trajectory, consistent 
with the persistence of slow-roll and the validity of the 
effective single field description [TJ1 [J21 [T51 QJj] . In this 
way, we offer a uniquely discriminable signature of the 
effect of higher dimensional operators that become rele- 
vant during inflation, the positive observation of which 
would allow us to infer properties of the parent theory in 
which inflation is embedded. 
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